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Abstract 

A gluon resonance G of mass below 1 TeV could be the origin of the tt forward- 
backward asymmetry observed at the Tevatron provided that new decay modes G — )• 
qQ, with q a standard quark and Q its massive excitation, make G broad enough. We 
consider all the different cases, with q the top, the bottom or a light quark and dominant 
decay modes Q — )• Wq' or Q — t- Zq. We show that current experimental searches are 
unable to probe the model, but that minimal departures from these analyses can explore 
a large region of its parameter space for the current LHC luminosity. This includes the 
challenging case with the new quarks decaying mostly into light quark flavors. In some 
channels not only the heavy quark but also the massive gluon can be reconstructed, 
which would stablish the origin of the ti asymmetry. Similar analyses can be applied 
to more general models with new massive gluons and vector like quarks. 



1 Introduction 



With over 10 fb~^ of recorded data at the Tevatron and more than 4 fb~^ at the Large Hadron 
CoUider (LHC), physics beyond the standard model (SM) is currently being searched with 
a very important degree of detail. Until now no discovery has been reported by any exper- 
imental collaboration, and bounds on many extensions of the SM rise up to the TeV scale 
and sometimes higher. As an alternative, these results may just imply that the experimental 
signature of the new physics is peculiar and easy to miss despite being at relatively low 
scales, as preferred by naturalness arguments. In this article we take this approach to study 
the tt forward-backward (FB) asymmetry at the Tevatron [Ij, arguably the most intriguing 
departure from the SM predictions. We show that it can be explained by new physics below 1 
TeV that could be difficult to see in the first round of experiments unless dedicated analyses 
are considered. 

Due to the relatively large coupling of the top quark to the electroweak (EW) symmetry 
breaking sector, new physics stabilizing the latter could also appear in top-quark observables. 
This generic argument makes the 2 — 3 a deviation versus the standard value in the Tevatron 
asymmetry specially interesting. Even if it is not statistically significant at the level of 
discovery, the consistency among different CDF and DO measurements strengthens the case 
for new physics. However, any candidate responsible for the asymmetry has to be carefully 
disguised, as its large contribution there should not translate into any significant departure 
from the SM in other related observables. In particular, the tt total cross section, its invariant 
mass distribution, dijet production, same sign top production, or the tt charge asymmetry 
at the LHC are observables where correlated anomalies could be expected |2]. 

In a recent work [3j we have shown that an s-channel gluon resonance G of relatively 
low mass [Mq ^ 1 TeV) could explain the large value of the asymmetry consistently with 
all the other observations. (See for alternative explanations of the Tevatron asymmetry 
in terms of massive gluons and ways to discover them.) It should have small-close to axial 
couplings to the light quarks {g1 ~ — fl'^) together with a large coupling to the right-handed 
top quark, features that are obtained in Higgsless models after imposing consistency with 
EW precission data [5]. The key ingredient would be a large gluon width, Tq = (0.5- 
0.7)Mg, provided by new decay modes of type G — > Qq, qQ, where g is a standard quark 
and Q a massive vectorlike excitation. In composite holographic models these fields can be 
understood as Kaluza-Klein modes of the standard quarks. The large gluon width in this 
framework requires a proper treatment of energy- dependent effects. In particular, a Breit- 
Wigner with constant width would offer a poor description of the gluon-mediated amplitude. 
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where gy'^ = {gn ± (7l)/2 are the vector and axial coupling of the massive gluon to Q and g, 
respectively. The large width will then reduce all gluon effects at -\/I > mq + (like a peak 
in the tt or the dijet distributions) while leaving unchanged lower energy effects (namely, 
the Tevatron FB asymmetry). 

For this mechanism to properly explain the asymmetry we need 700 GeV < Mq < 
900 GeV and 400 GeV < mq < 700 GeV [3]. The low masses of the gluon and the new 
quarks, together with the sizable couplings required to generate the large width, make the 
production of single new quarks mediated by the massive gluon a very attractive channel 
at the LHC. In this article we investigate its potential to probe this scenario. The signal 
there will depend strongly on the nature of the vectorlike quark involved. In section [2] we 
classify all the possibilities and introduce a benchmark model that provides contributions in 
all the different channels. In section |3] we study single vectorlike quark production involving 
the third generation, and in section |4] we discuss the channels with light flavors. In both 
cases we show that current analysis could easily miss the model, whereas specific searches 
would very likely reveal the mechanism responsible for the Tevatron asymmetry. Section [5] 
is devoted to our conclusions. 



2 A benchmark model 

In this section we introduce a benchmark model that successfully reproduces the Tevatron 
FB asymmetry with no conflict with other experimental tests. It contains simultaneously 
all possible decay channels and, therefore, allows us to perform a comprehensive study of 
the stealth gluon scenario. The model admits variations where one or several channels are 
suppressed while the others are enhanced in such a way that the total gluon width does not 
change significantly. We take Mq = 850 GeV, although similar setups can be obtained for 
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gluon masses as low as 700 GeV. We fix the couplings to G of the SM quarks to 

gl = 0.3 5„ 5^ = <7^ = -0.3^„ g'R = +Ags, g\ = g\ = ^, (3) 

where gg is the strong coupling constant. For the vectorlike quarks, we will assume the 
presence of six fields, corresponding to the excitations of tji^ bu and the four hght flavors ql 
(from now on we use q and Q just for the {u, d, s, c) quarks and their excitations). We fix 
their masses to 

Mt = 450 GeV, Ms = Mq = 600 GeV, (4) 
and their flavor-changing couplings to the heavy gluon to 

9r =4^,, gi'^3.5gs, g'i' = 3.5 g^ . (5) 

These values imply a total width Vq ~ 0.7 Mq and the decay branching fractions 

BR{G ^tt)^ 0.2, BR{G Tt, tf) ^ 0.24, 

BR(G ^ Bb, bB) ^ 0.11, BR{G Qq, qQ) ^ 0.44. (6) 



As we mentioned above, the benchmark model just defined has the advantage that all 
possible channels are present simultaneously. However, when studing the possible collider 
implications of this scenario we will also consider the extreme cases where all but one G 
decay modes are absent: 



Extreme T model: 
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and all the other couplings unchanged. In these cases the heavy gluon has a 20% branching 
ratio into tt and 80% into the new channel. Note that in some of these models the required 
coupling is unrealistically large. We just take them as limiting examples to get clear idea 
of the LHC reach for these signatures (reahstic models should he somewhere in between the 
benchmark and the extreme cases). 

The new heavy quarks will then be produced through G in the s-channel as QQ pairs 
or as a single particle together with a standard quark, Qq. Pair production will also receive 
the standard QCD contribution (in fact, due to the axial nature of the G coupling to light 
quarks, the interference terms cancel and away from the resonance pair production is like 
in the SM). Single heavy-quark production, on the other hand, is unsuppressed and opens 
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Figure 1: rriu distribution at the Tevatron for 5.3 fb~^ in the SM (sohd blue), the benchmark 
model (points with error bars) and the extreme T case (dotted black). We include the 
contribution from Tt, tT and Bb, hB when present. 



kinematically at lower energies (a/S = niq + mq <^ Smg), appearing as a very promising 
mechanism unexplored in previous literature. The vectorlike quarks will then decay in a 
model-dependent way, according to their electroweak quantum numbers and their mixing 
with the SM quarks. Assuming weak couplings, their width will be narrow, and a simple 
scaling allows to go from one model to another. To be definite we will take the branching 
ratios obtained in the large-mass limit of the usual Higgsless models, 

BR(Q ^ W^g') = ^ , BR(g^Zg) = ^. (10) 

Higgs decays can potentially lead to interesting signatures |i.6] but we defer the corresponding 
analysis to future work. In this paragraph, we have denoted with Q the six vectorlike quarks. 

With these assumptions the final states produced in qq collisions will be the following 
(the conjugated processes are not explicitly shown but are included in our analyses): 

(i) W+W'hh, from 

qq ^ G ^ Tt ^ {W+h)W-h (11) 

and 

qq G ^ Bb ^ {W-t)b ^ {W-W+b)b . (12) 

(a) Ztt , from 

qq^G^Tt^ {Zt)t {ZW^b)W-b . (13) 
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( in ) Zbb , from 



qq-^G^Bb^ {Zb)b. 



(14) 



(iv) W+iets, from 



qq ^ G ^ Qq ^ {Wq)q 



(15) 



(v) Z+jets , from 



qq ^ G ^ Qq {Zq)q . 



(16) 



In the next two sections we show that these signals do not introduce observable anoma- 
lies in current LHC analyses, but that simple modifications in the reconstruction of the final 
state could very likely provide a signal. The impact of this scenario on top-quark physics at 
the Tevatron has been discussed in [3], where we name it as the stealth gluon model due to 
its ability to explain the FB asymmetry without introducing anomalies (peaks or tails) in 
the ti invariant mass distribution [rritt)- In particular, it implies A^^{mtt < 450 GeV) = 0.12 
and A^^{mu > 450 GeV) = 0.33, values that are compatible with the DO and CDF obser- 
vations [Ij. The rritt distribution at the Tevatron is given in Fig. [H where we compare the 
reconstruction as tt pairs of all the events giving W^W~bb in the benchmark model with 
the SM prediction. In our simulation we have followed the analysis in [7] and have gener- 
ated the events with MADGRAPH/MADEVENT v4 |8] (with the matrix element properly 
modified to include the energy dependence of the width), using PYTHIA [10] for hadroniza- 
tion/showering effects and PGS4 [11] and DELPHES 1.9 [12] for detector simulation. We 
include in the figure the prediction in the extreme T model (the prediction in the extreme 
B model is similar, whereas in the extreme Q model it is below the benchmark one). The 
deviations are never larger than 2.5 cr (assuming statistical errors only), and below 2 a in 
all the bins for the benchmark and the extreme-Q models. 

3 Single T and B quark production at the LHC 
3.1 W+W-hh channel 

As described in the previous section, the new processes qq — ?■ Tt, Bb followed by the charged- 
current decay of the heavy quark will result in the same W~^W~bb final state as ti production. 
In our model this signal would add to the one from top-quark pairs produced through the 
massive gluon, and it is then necessary to check that these processes do not imply any 
observable excess in current analyses of tt production or fourth generation TT searches. 
Measurements of the tt mass invariant distribution at the LHC have been reported in [13] . 
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We have simulated the analyses in the first two works of this reference and studied the effect 
of the channels 

pp Tt, tf, Bb, bB (17) 

together with all the contributions to tt production. We show the result in Fig. [2] (we have 
followed the second Ref. in [13] for 0.2 fb^^; the third reference uses the dilepton channel 
and a larger data set, implying a very similar sensitivity). In the plot we have assumed a 
10% uncertainty in the ti prediction and allowed a normalization factor (within this 10%) 
to correctly reproduce the three bins around the peak at rriu ~ 500 GeV. We show the SM, 
the benchmark model (with statistical error bars) and the extreme T model. The deviation 
in the extreme B case is similar to the one in the extreme T model, whereas the extreme Q 
case is closer than the benchmark to the SM. The ^ 20% excess at rriu = 600-900 GeV in 
the extreme T and B models seems in the limit of being probed with the current LHC data. 
Increasing the luminosity to 4 fb~^ we find 8 consecutive bins with differences above 3 a for 
the DELPHES simulation and 7 consecutive ones for the PCS simulation in the case of the 
extreme T model. The benchmark and extreme Q models are not as clear. For instance, 
using PGS we find 3 and 2 consecutive bins with departures larger than 3 a in these cases for 
a luminosity of 4 fb~^ (in all our estimates we only include statistical errors). In summary, 
in our model one could expect a 10% excess relative to the SM prediction in all the bins 
below 1 TeV. These events are just ti pairs mediated by the heavy gluon G. In addition, the 
bins between 600-900 GeV could be increased an extra 15% with Tt and/or Bb events that 
are reconstructed as ti pairs. 

Another LHC study sensitive to our model is the search for a fourth generation of quarks 
produced as TT pairs [H]. We have reproduced the corresponding CMS analysis and plot 
our results in Fig. |3]for the muon channel with the published luminosity of 0.821 fb~^. Our 
results are similar to the ones obtained for ti production. The benchmark and the extreme 
Q models are not visible, whereas the extreme T and B models are starting to be probed 
by the data. We plot in Fig. 3 the SM, the benchmark and the extreme T cases in solid 
blue, data points (with error bars), and dotted black, respectively. The left panel shows 
the Ht distribution (defined in this case as the scalar sum of the pr of the jets, the charged 
lepton and the missing Et), and the right panel gives the T reconstructed mass in the events 
generated with our model(s) and with the SM. In both plots the number of standard events 
has been normalized by the same factor. We have also checked that pair production of T 
quarks give in our model a negligible contribution, compatible with the bound obtained 
in [ll]. Similarly, the recent search for pair production of vectorlike T quarks decaying to 
Zt [15] does not imply any restriction to our model. 

Our results indicate that the model, proposed to explain the large FB Tevatron asymme- 
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Figure 2: rriu distribution at the LHC for 0.2 fb^^ in the SM (sohd blue), the benchmark 
model (points with error bars) and the extreme T model (dotted black). We include the 
contribution from Ti, tT and Bb, bB when present. 

try, is almost invisible in ti — )■ W~^W~bb searches. The reason for that is twofold. First, the 
large gluon width suppresses the number of ti events in the region mu = 600-900 GeV, while 
its axial couplings to the light quarks does the same at lower and higher invariant masses. 
Second, Ti or Bb events are reconstructed as ti or TT pairs, resulting into a poorer fit and 
a wider spread. The key to isolate events of type Ti would be to reconstruct them not like 
two objects with the same mass, but like a t quark plus a T quark of arbitrary mass. These 
events will only occur at large invariant masses, rriTt > + mt, a. region already accessible 
at the LHC with the current luminosity. Therefore, we can use the more stringent cuts used 
in the TT analysis of [IA\ (we use the muon channel). Actually, we will require the hardest 
jet to have px > 200 GeV instead of the 120 GeV of that reference. We will then identify 
just one 173 GeV t quark (using a similar to the one used in the first reference of [13] and 
requiring < 10) and will plot the mass of the second one in events of invariant mass above 
600 GeV (Fig. HJ left panel) for SM and extreme T model simulations. We have normalized 
the plots to the recorded luminosity of 4 fb~^. As it is apparent in the plot, we find three 
consecutive bins around rriT = 450 GeV departing more than three sigmas from the SM 
prediction even in the benchmark model. Counting the total excess S of events versus the 
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Figure 3: TT search at the LHC for 0.821 fb~^. Left paneh Ht distribution. Right paneh 
rrifit distribution. In both cases we show the predictions in the SM (sohd blue), in the 
benchmark model (data points with statistical errors) and in the extreme T case (dotted 
black). We include the contribution from Tt,tT and Bb,bB when present. 



standard background B on the peak (three bins between 350 and 500 GeV) we get 

5* J 8, benchmark, 

. — ~ < (18) 
I 21, extreme T. 

Thus, the extreme T case would imply a stunning deviation in this kind of searches, and 
even the benchmark model could show evidence for new physics. With the large excess in 
the extreme T model one can also try to reconstruct the massive gluon peak. In order to do 
that, we remove the total invariant mass cut and compute the total invariant mass rriTt for 
the events with a reconstructed T mass above 350 GeV. The result is shown in Fig. HJ right 
panel. Although the SM and the new physics model peak in the same region, the factor of 
~ 3(2) excess in the extreme B (benchmark) model is quite evident. 

The Bb, bB — W^W~bb channel is slightly different. Instead of producing two top-like 
objects, the heavy bottom decays into a W plus a top that subsequently decays into another 
W (with opposite charge) and a b. We will still follow the selection procedure in our previous 
analysis, with the cuts in [2] (muon channel) except for the cut on the pt of the hardest jet, 
that is moved from 120 GeV to 200 GeV and a < 10 (again we use a similar to the one 
used in the first reference of [13]) choosing the best configuration reconstructing a 173 GeV 
top quark, and will plot the invariant mass of this t quark plus the extra W. The result is 
shown in Fig. [5] for the benchmark and the extreme B models with two different cuts in the 
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Figure 4: Left panel: Reconstruction of at the LHC Right paneh Reconstruction of ttlg- 
In both cases we have normahzed the distributions to 4 fb~^ data and represent the results 
for the SM (solid blue), the benchmark model (data points with statistical errors) and the 
extreme T case (dotted black). Details of the reconstruction method can be found in the 
text. 

total invariant mass distribution and 4 fb~^. In this case, our reconstruction of the B quark 
is not as clear as the one of the T quark, and more sophisticated analyses should be used to 
dig out the signal from the background. Nevertheless, we will see in the next section that 
the extreme B model can be probed much more efficiently using the neutral decay of the B 
quark. 



3.2 Zhh channel 

Let us now turn to the neutral decays of the heavy T and B quarks, starting with the 
Bh, bB channel into a Zbb final state. The SM irreducible background to this process is 
small {cr{Zbb) with a leptonic Z decay is around 2 pb), whereas the background from final 
states with larger cross sections like Z+jets and ti can be reduced with a very simple set 
of cutslll To isolate the signal we will require two same-fiavor opposite-sign leptons with 
Pt > 25 GeV and |m/+;~ — mz\ < 25 GeV, and two 6-tagged jets with pt > 20 GeV and 
|?7| < 2.8. We will also impose a veto on missing energy Ej- < 40 GeV, to reduce the tt 

^We have also checked that our model does not conflict with current searches oi H ^ ZZ Zhh [IH] or 
measurements oi Z + h cross-section [17] . 
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Figure 5: Reconstruction of tub at the LHC for 4 fb^^ in the SM (sohd blue), the benchmark 
model (data points with errors) and the B case (dotted black). We consider the cuts rriBb > 
600 GeV (left) and m^j > 700 GeV (right). Details of the reconstruction method can be 
found in the text. 

background. With this selection we compute the invariant mass of the Z and the hardest of 
the two 6-jets (denoted by bh), since the b quark from the decay of the heavy B is typically 
the hardest one. We plot the result in Fig. |6l In the left panel we show the mzb^ invariant 
mass distribution in the SM, the benchmark model and the extreme B case. It is clear 
that the distributions in the SM and the new model peak in very different regions. The 
benchmark model leads to too small a cross section and would require higher luminosity 
for discovery. The extreme B model, however, shows a clear peak with a total number of 
40 events at mzbh ~ = 600 GeV, versus ~ 3 background events, implying a statistical 
significance of 

S 

21, {Zbb for extreme B). (19) 



^B 

Given the presence of a distinct peak we can attempt to reconstruct the mass of the heavy 
gluon. In the right panel of Fig. |6]we show the total invariant mass of the three objects Zbb 
for the events passing the cuts. Due to the large width of the heavy gluon (the kinematical 
threshold prevents the full width to be apparent at energies below ~ 600 GeV) the number 
of events peaks slightly below Mg = 850 GeV, but the effect is clearly observable. The 
approximate statistical significance of the excess above 600 GeV is 

S 38 

—= = 17, {Mg peak in Zbb for extreme B). (20) 
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Figure 6: Left panel: reconstruction of mzbh at the LHC Right paneh reconstruction of 
m^i,l to show the heavy gluon mass. In both cases we have normahzed the distributions to 4 
fb~^ of data and have represented the SM with thick sohd blue line, the benchmark model 
with thin solid red line and the extreme B case (data points with statistical errors). 

The Zbb channel appears then as very promising even with the very simple cuts that 
we have used. In the extreme case the reconstruction of the B quark and of the massive 
gluon at the 4 fb~^ LHC could be correlated with the ti anomalies discussed in Section 3.1, 
disentangling the origin of the Tevatron FB asymmetry. 

3.3 Zti channel 

The Ztt production channel resulting into a ZW~^W~bb final state has also a very small SM 
background, but it is harder to reconstruct due to its large multiplicity. Instead of trying to 
reconstruct the T mass, it is simpler to reconstruct the total final state in the search for the 
massive gluon. We do that requiring (i) three charged leptons with px > 25 GeV, and at 
least two of them with the same flavor and opposite sign reconstructing the Z within 25 GeV; 
(ii) at least two 6-tagged and at least two non-6-tagged jets with px > 20 and \ri\ < 2.8. 
We reconstruct the neutrino momentum using the on-shellness condition for a W and take 
the two hardest jets and 6-jets if there are more of them. The result is shown in Fig. [71 The 
extreme T model shows a clear peak with ^ 36 events with no expected background events 
(the benchmark gives a weaker deviation). A more detailed analysis, trying to reconstruct 
both top quarks, would certainly help in the reconstruction of the heavy T mass. Since the 
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Figure 7: Total invariant mass reconstruction for the Zti channel in the SM (solid blue), 
benchmark (solid red) and extreme T (data with statistical errors shown as a band) models 
for the Ztt analysis described in the text for the LHC with 4 fb~^. 

extreme T model would also show up in the charged decay channel, a hint on the T mass 
could be used in the reconstruction of this channel. 

4 Light flavor excitations: Wq'q and Zqq 

We have seen in previous sections that the production of single T or 5 quarks tend to 
introduce anomalies in current searches and could be seen if the reconstruction algorithms 
are slightly modified. However, Qq production is less apparent in these searches, being the 
best example of stealth new physics [3J. We discuss in this section the best strategy to observe 
the extreme Q model at the LHC. In the benchmark (extreme Q) model the production of 
heavy excitations Q of the light flavors has a total cross section of 2.9 (5.4) pb at the 7 
TeV LHC, resulting with a 2:1 ratio the final states Wq'q and Zqq. The SM irreducible 
background is 17 nb for W plus > 1 jets and 6 nb for Z plus > 1 jets. Therefore, we 
need to impose stringent cuts to disentangle our signal from these large backgrounds. First 
of all, these extra Qq events will only appear at invariant masses above rriQ = 600 GeV, 
with the maximum at ^ 700 GeV. In addition, the jet from the decay of the heavy quark, 
with a ~ will be typically harder than the second jet. We should then impose 

an stringent cut on the hardest jet in order to reduce the SM backgrounds. In particular. 
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Figure 8: Left panel: transverse mass for the Wjh system in the Wjj analysis described in 
the text for the SM (solid blue), benchmark model (data points with errors) and extreme 
Q model (dotted black). Right panel: Result of the fit of the mzj^ distribution for the Zjj 
analysis described in the text for the SM (solid blue), extreme Q model (data points with 
statistical errors) and the fit to both distributions (dotted black). Both plots are for the 7 
TeV LHC with 4 fb'^ 

requiring a hardest jet with px > 150 GeV on top of the cuts defined in Ref. [18] reduces 
the W+jets background to manageable levels. We show in Fig. [S] (left panel) the transverse 
mass distribution of the W and the hardest jet. The signal does not seem significant in the 
benchmark model but may be observable in the extreme Q case, with 6 bins departing more 
than 3 standard deviations from the expected background. 

The neutral case is even more promising. Requiring two same-fiavor, opposite-charge 
leptons with px > 25 GeV that reconstruct the Z mass within 25 GeV, and two or more 
jets, with pt > 150 GeV for one of them, and computing the invariant mass of the Z and 
the hardest jet, we obtain the distribution in Fig. |8] (right panel). Although the benchmark 
model is still unobservable, there is a clear peak for the extreme model. We have fitted the 
signal plus background histogram to a Crystal Ball plus gaussian shape and obtained an 
excess of 170 events over the expected 540 background events in the region of two stardard 
deviations around the center of the gaussian. This leads to a statistical significance of 7a 
and a best fit of itlq = 590 GeV, very close to the actual heavy quark mass. This analysis is 
interesting as it gives a very clean signal for a model that is otherwise very difficult to find. 
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5 Summary and discussion 



The Tevatron A** anomaly is a strong motivation for a search of correlated effects from new 
physics at the 7 TeV LHC An explanation with physics above ~ 1 TeV seems disfavoured 
by data p]. If it is below this energy, then a large width could be the key reason why it 
has escaped detection in the usual observables. In a recent work we proposed that a massive 
gluon with new decay modes G — )■ Qq would be a promising candidate, and here we have 
studied in some detail the consequences of these processes in current analyses and possible 
signals to be searched. 

ATLAS and CMS are studying tt and TT production. We have shown that Tt production 
could also be explored just by slightly changing the criteria of reconstruction. The channel 
Bb provides the same WWbb signal and could also be studied there. 

We have also discussed new channels that, if analyzed, could reveal single heavy quark 
production at the LHC. In particular, Zqq where the Q quark is reconstructed with the 
Z boson and the highest-energy jet looks promising. Other signals, like Zbb or Ztt, are 
predicted here and have small SM backgrounds. The study of these channels is well motivated 
by holographic models. If the Tevatron anomaly is due to new quark interactions below the 
TeV, then they should be searched for at the LHC, since there seems to be few hideouts 
for the new physics beneath. We have focused our study on the region motivated by the 
Tevatron asymmetry but our analyses can be also applied to a wider range of couplings and 
quark and gluon masses. 

Note added 

During completion of this work we became aware of Refs. [19] in which similar ideas to the 
ones presented here are being investigated. See also [20\ . 
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